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What do dynamical cluster masses really tell us about dynamics? 
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The diagnostic age versus mass-to-light ratio diagram is often used in attempts to constrain the shape of the stellar 
initial mass function, and the stability and the potential longevity of extragalactic young to intermediate-age massive star 
clusters. Here, we explore the pitfalls associated with this approach and its potential for use with Galactic open clusters. 
We conclude that for an open cluster to survive for any significant fraction of a Hubble time (in the absence of substantial 
external perturbations), it is a necessary but not a sufficient condition to be located close to the predicted photometric 
evolutionary sequences for 'normal' simple stellar populations. 

© 2008 WILEY- VCH Vcrlag GmbH & Co. KGaA, Wcinhcim 



1 A diagnostic diagram? 

Over the past few years, detailed studies of the stellar con- 
tent and longevity of extragalactic massive star clusters have 
increasingly resorted to the use of the age versus mass-to- 
light (A// L) ratio diagram as a diagnostic tool (see Fig.[T]for 
an up-to-date version), where one usually compares dynam- 
ically determined M/L ratios with those predicted by the 
evolution of 'simple' stellar populations (SSPs; e.g.. Smith 
& Gallagher 2001 ; Mengel et al. 2002; McCrady et al. 2003, 
2005; Larsen et al. 2004; Bastian et al. 2006; Goodwin & 
Bastian 2006; Moll et al. 2008; see de Grijs & Parmentier 
2007 for a review). Based on high-resolution spectroscopy 
to obtain the objects' integrated velocity dispersions, a, and 
on high spatial resolution imaging to obtain accurate pro- 
jected half-light radii, rhi, most authors then apply Spitzer's 
(1987) equation, 
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to obtain the dynamical cluster masses, Mdyn (G is the grav- 
itational constant and 77 w 9.75 is a dimensionless param- 
eter which is usually assumed to be constant; but see Fleck 
et al. 2006; Kouwenhoven & de Grijs 2008). 

Despite a number of simplifying assumptions (see, e.g., 
de Grijs & Parmentier 2007 for a review; Moll et al. 2008), 
one can get at least an initial assessment as to whether a 
given cluster may be (i) significantly out of virial equilib- 
rium, in particular 'super- virial', (ii) significantly over- or 
underabundant in low-mass stars, or (iii) populated by a sig- 
nificant fraction of binary and higher-order multiple sys- 
tems. This has led to suggestions that, in the absence of 
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Fig. 1 Updated version of the M/L ratio versus age di- 
agnostic diagram for young and intermediate-age massive 
star clusters. The numbered data points follow de Grijs & 
Parmentier (2007; their fig. 2); overplotted are the Maraston 
(2005) SSP predictions for a Salpeter (1955) and a Kroupa 
(2001) stellar IMF, as well as the effects of varying effec- 
tive star-formation efficiencies, based on Goodwin & Bas- 
tian (2006). 



e-mail: R.deGrijs@sheffield.ac.uk 



®W]LEY 

JiTiterScieTice* 



© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim 



2 



R. de Grijs et al.: Dynamical cluster masses 



significant external perturbations, young massive clusters 
(YMCs) located in the vicinity of the SSP models and aged 
^ 10^ yr may survive for a Hubble time and eventually be- 
come old globular cluster (GC)-like objects (e.g., Larsen et 
al. 2004; Bastian et al. 2006; de Grijs & Parmentier 2007). 

In this contribution we will first highlight the uncertain- 
ties associated with a few of the YMCs (Section 2) - and 
thus question the use of the velocity dispersion as a tracer 
of their gravitational potential - and then proceed to apply 
the same diagnostic approach to a small sample of Galac- 
tic open clusters (Section 3). In Section 4 we will provide a 
summary, emphasising the pitfalls associated with the use of 
dynamical mass estimates as a proxy for cluster longevity. 

2 Closer inspection 

Let us now examine a few of the YMCs from Fig.[T]in more 
detail. Using the dynamical mass estimate from Mengel & 
Tacconi-Garman (2007), combined with the integrated pho- 
tometry of Piatti et al. (1998), in de Grijs & Parmentier 
(2007) we concluded that the Galactic YMC Westerlund 1 
appears to be consistent with a normal Kroupa (2001) or 
Salpeter (1955)-type initial mass function (IMF), despite 
the significant uncertainties in the observables. Since in the 
V band, on which the Piatti et al. (1998) photometry was 
based, the confusion between the cluster members and the 
Galactic field stellar population is substantial (because of 
the significant extinction along this sightline), we obtained 
imaging observations in the / band, where this confusion is 
significantly reduced (de Grijs et al. 2008). 

The combined integrated magnitude of the three bright- 
est red supergiants, yellow hypergiants, and blue supergi- 
ants is some 40% of the cluster's total integrated /-band 
flux (see de Grijs et al. 2008 for details). These nine sources 
stand out from the overall stellar luminosity function, which 
appears to otherwise have been drawn from a 'normal' IMF. 
Each of these sources is in a rare, short-lived phase and so 
the luminosity of the cluster might be expected to vary sig- 
nificantly on short time-scales. Thus, this serves as a clear 
caution that stochasticity in the cluster's IMF (e.g., Brocato 
et al. 2000; Barker et al. 2008), as well as stochasticity in the 
numbers of stars in unusually luminous post-main-sequence 
evolutionary stages (e.g., Cervino & Valls-Gabaud 2003; 
Cervinno & Luridiana 2006) may contribute significantly 
to variations in a cluster's AI/L ratio. 

Thanks to Bastian & Goodwin (2006) and Goodwin & 
Bastian (2006), we are now confident that we understand 
why some of the YMCs deviate significantly from the ex- 
pected evolution - for a 'standard' IMF - of the stellar M/L 
ratio at ages of up to a few x 10^ yr However, the location 
of the massive cluster F in M82 remains to be understood, 
despite having been the subject of a considerable number of 
recent studies (see, e.g., Bastian et al. 2007 for an overview). 

Smith & Gallagher (2001) concluded that, assuming that 
their age (60 ±20 Myr), mass (1.2 ±0.1 x IO^Mq) and lu- 
minosity (My = -14.5 ± 0.3 mag; Lv/M^yn = 45 ± 



13 Lv,q/Mq) measurements for the object were correct, it 
must either have a low-mass cut-off of its mass function at 
^ 2 — 3Afo, or a shallow mass function slope, a ^ 2 for a 
mass function including stellar masses down to 0.\Mq. Mc- 
Crady et al. (2003) suggested, based on near-infrared pop- 
ulation synthesis modeling, that Smith & Gallagher (2001) 
had overestimated their ages as well as the half-mass radius 
of the cluster; they derived the latter more accurately based 
on higher-resolution imaging. However, in a recent paper 
Bastian et al. (2007) used optical spectroscopy, confirming 
the YMCs age to be in the range of 50-70 Myr. However, 
contrary to the earUer studies, Bastian et al. (2007) conclude 
that M82-F is subject to a large amount of differential ex- 
tinction, thus rendering earlier luminosity estimates based 
on the assumption of foreground-screen extinction very un- 
certain. They conclusively show that the apparently large 
degree of mass segregation in the cluster, derived by Mc- 
Crady et al. (2005) based on the variation of 7-hi as a func- 
tion of wavelength, is in fact caused by this differential ex- 
tinction. Applying their spatially resolved extinction correc- 
tions to Hubble Space Telescope imaging data, Bastian et al. 
(2007) also suggest that the original cluster size reported by 
Smith & Gallagher (2001) was underestimated. Using their 
updated (larger) size estimate, the discrepant M82-F data 
point moves closer to the SSP predictions, hence suggest- 
ing that the cluster may be less enigmatic than previously 
thought. 



Finally, we single out YMC NGC 1 140-1 (Hunter et al. 
1994; de Grijs et al. 2004; Mofl et al. 2008) for a closer 
look. There are various potential explanations for the differ- 
ence between its observed and the SSP model M/Ly ratio. 
Assuming that this difference is due solely to the cluster 
being out of virial equilibrium, its effective star-formation 
efficiency (eSFE; Goodwin & Bastian 2006) is around 10- 
20%. The eSFE is a measure of the extent to which a clus- 
ter is out of equilibrium after gas expulsion, on the basis 
that the virial ratio immediately before gas expulsion was 
Qvii = T/\n\ = l/2(eSFE), where T and Vl are the ki- 
netic and potential energy of the stars, respectively, and a 
system in virial equilibrium has Qvir — 1/2. The eSFE 
corresponds to the true SEE if the stars and gas were ini- 
tially in virial equilibrium. Goodwin & Bastian (2006) pre- 
dict that a cluster with an eSEE of 10% will lose 80% of its 
mass in 20 Myr and not remain bound. As such, the cluster 
would disperse over a relatively short time-scale and, there- 
fore, not evolve into a second-generation GC. More likely, 
however, the crowded nature of the region may have caused 
the cluster velocity dispersion (and hence the M/L ratio) to 
be overestimated. The broadening of the red-supergiant fea- 
tures seen in the knot spectra may not have been due solely 
to the virial motions of the brightest cluster within the knot, 
as was assumed, but may have had contributions from the 
motions of nearby contaminating clusters instead (cf. Moll 
et al. 2008). 
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Fig. 2 Westerlund 1 in the diagnostic age versus MjLy 
ratio diagram. The open circle represents the shift in the 
cluster's locus from its current position (solid bullet) if we 
were to exclude the nine brightest stars; this exemplifies the 
uncertainties introduced by stochastic IMF sampling and by 
having fortuitously caught the cluster at a time when it is 
dominated by a small number of very bright stars. 




6 7 8 9 10 

log( Age / yr ) 

Fig. 3 As Fig. 1, but for those Galactic open clusters 
for which velocity dispersion measurements are available. 
Numbered clusters: 1, NGC 1976 (Orion Nebula Cluster); 
2, NGC 2168; 3, NGC 2516; 4, NGC 2632; 5, NGC 2682; 
6, NGC 3532; 7, NGC 5662; 8, NGC 6705; 9, Pleiades; 10, 
Coma Berenices; 11, Hyades. 



3 Galactic open clusters 

We will now explore whether we can also use the same diag- 
nostic diagram to assess the stability, formation conditions, 
binarity and longevity of those open clusters in the Milky 
Way for which the requisite observational data exist in the 
literature. Using the observational data compiled in de Grijs 
et al. (2008), we appHed Spitzer's (1987) equation, Eq. O, 
to derive the dynamical masses for each of our sample clus- 
ters, and then calculated their MjLy ratios. The clusters' 
loci in the diagnostic diagram are shown in Fig. [3] 

We have also included the expected evolution of clus- 
ters formed with a variety of eSFEs. Owing to the nature of 
our sample, onl}Q the Orion Nebula Cluster (ONC; cluster 
1 in Fig. O is young enough so as to possibly be affected 
by the effects of rapid gas expulsion, as shown by the extent 
(in terms of age) of the long-dashed lines in Fig. [3] The ma- 
jority of our sample clusters are old enough (J^ 40 Myr) to 
have re-virialised after gas expulsion. The dynamical state 
of these objects is therefore dominated by the combined ef- 
fects of (internal) two-body relaxation, binary motions, and 
external perturbations. 

We note that, despite the sometimes significant uncer- 
tainties (parameter ranges) associated with the individual 

' Despite the extent of the error bar associated with the age estimates of 
the Pleiades, Kroupa et al. (2001) showed this cluster to have re-virialised 
by an age of 50 Myr, so that it is unlikely affected by the aftermath of the 
gas-expulsion phase. 



quantities required to calculate the dynamical M / Ly ratios, 
the sample clusters follow the general trend predicted by the 
SSP models rather closely. The fact that these clusters lie 
close to (although systematically somewhat below) the SSP 
predictions should not be a suprise. Clusters significantly 
below the SSP lines are dynamically 'hot' and expected to 
dissolve rapidly, whilst clusters significantly above the lines 
will be dynamically 'cold' and should (re-)virialise over a 
few crossing times to move closer to the canonical SSP 
lines. 

Most of the sample clusters are found somewhat below 
the SSP model curves. Errors in the core radii are expected 
to be random, and unbiased by the mass of a cluster. How- 
ever, the use of the core velocity dispersions and radii may 
introduce a systematic bias in the dynamical mass estimates. 
The majority of the star clusters in the sample are older than 
^ 10^ yr, which implies that they have ages greater than 
their (initial) half-mass relaxation times (cf. Danilov & Se- 
leznev 1994). Therefore, these clusters are expected to be 
close to energy equipartition, and thus are significantly mass 
segregated. Equipartition reduces the global velocity disper- 
sion of high-mass stars relative to low-mass stars, causing 
high-mass stars to migrate to the cluster core. Therefore, 
we might expect the core velocity dispersion of low-mass 
cores to underestimate the dynamical mass and thus pro- 
duce colder clusters - as observed. 

Alternatively, the effect of binaries within clusters may 
well account for most of the displacement of the observed 
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cluster positions below the model curves. Kouwenhoven & 
de Grijs (2008) pointed out that if the velocity dispersion of 
binary systems was similar to the velocity dispersion of the 
cluster as a whole, the observationally measured velocity 
dispersion would overestimate the mass of a cluster Based 
on a comparison with Kouwenhoven & de Grijs (2008), it 
appears that the vast majority of the sample clusters are in- 
deed expected to be binary dominated (see de Grijs et al. 
2008 for details). 

We also note that the cluster masses may well have been 
overestimated by factors of a few through the universal use 
of Eq. ([T]i. In particular, for highly mass-segregated clus- 
ters containing significant binary fractions, a range of stel- 
lar IMF representations, and for combinations of character- 
istic relaxation time-scales and cluster half-mass radii, the 
adoption of a single scaling factor 77 « 9.75 introduces sys- 
tematic offsets. To correct for these, we would need to adopt 
smaller values of 77 (e.g.. Fleck et al. 2006; Kouwenhoven & 
de Grijs 2008), and this would thus lead to dynamical mass 
overestimates if 77 = 9.75 were assumed. 

In addition, the mass functions (MFs) of clusters will be 
altered due to the preferential loss of low-mass stars dur- 
ing two-body relaxation. This will result in an increasingly 
'top-heavy' MF in clusters, which will raise them to lower 
M/Lv ratios than would be expected from the canonical 
SSP models. The degree to which the MF will change de- 
pends on the two-body relaxation time which, to first order, 
depends on the mass of the cluster (and also on its size; how- 
ever, we ignore this for now). Thus, low-mass clusters are 
expected to have top-heavy MFs as compared to high-mass 
clusters. Therefore, one would expect low-mass clusters to 
lie some way above the canonical SSP models, and high- 
mass clusters to lie slightly above these lines. However, the 
effect of binaries in clusters is to increase the observed ve- 
locity dispersion and so overestimate the masses of clus- 
ters, shifting them back down towards, and even below, the 
canonical SSP lines. 

4 Concluding remarks 

In this contribution we have explored the usefulness of the 
diagnostic age versus M / L-ratio diagram in the context of 
YMCs and Galactic open clusters. This diagram is often 
used in the field of extragalactic young to intermediate-age 
massive star clusters to constrain the shape of their stellar 
IMF, as well as their stability and the likelihood of their 
longevity. Using the massive young Galactic cluster Wester- 
lund 1 as a key example, we caution that stochasticity in the 
IMF introduces significant additional uncertainties. There- 
fore, the stability and long-term survival chances of Wester- 
lund 1 remain inconclusive. Similarly, we have highlighted 
observational issues associated with M82-F and NGC 1 140- 
1 which render the use of their velocity dispersions as a di- 
agnostic tool questionable at best. 

For low-mass open clusters it seems clear that the effect 
of binaries, mass segregation, and the dynamical alteration 



of mass functions by two-body relaxation are important fac- 
tors that cannot be ignored. Most importantly, however, we 
conclude that for an open cluster to survive for any signif- 
icant length of time (in the absence of substantial external 
perturbations), it is a necessary but not a sufficient condition 
to be located close to the predicted photometric evolution- 
ary sequences for 'normal' SSPs. This is based on two of 
the sample clusters, the Coma Berenices cluster (number 10 
in Fig. 3) and the Hyades (11), which are known to be in 
a late stage of dissolution, yet lie very close indeed to ei- 
ther of the evolutionary sequences defined by the Salpeter 
(1955) or Kroupa (2001) IMFs. 
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